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We calculate the self-energy of two-dimensional fermions that are coupled to transverse 
gauge fields, taking two-loop corrections into account. Given a bare gauge field propa- 
gator that diverges for small momentum transfers q as l/g* 7 , 1 < r\ < 2, the fermionic 
self-energy without vertex corrections vanishes for small frequencies u) as E(u) oc u) 1 
with 7 = < 1. We show that inclusion of the leading radiative correction to the 

fermion - gauge field vertex leads to T,(u>) oc lu"'[1 — a v ln(uiQ / u>)] , where a v is a pos- 
itive numerical constant and ujq is some finite energy scale. The negative logarithmic 
correction is consistent with the scenario that higher order vertex corrections push the 
exponent 7 to larger values. 

1. Introduction 

The problem of two-dimensional fermions that are coupled to a transverse gauge 
field has recently received a lot of attention, because it arises as effective low-energy 
theory in two different physical contexts. On the one hand, the infrared physics 

described by 
The other 



I 



of the two-dimensional t — J-model has been argued to be corre^ 
fermions and bosons that are coupled to an Abelian gauge field! 
example are half-filled quantum Hall systemsu, where a fictitious Chern-Simons 
gauge field can be used to attach two quanta of a magnetic flux to the physical 
electron, thus forming a new (presumably stable) quasi-particle, the so-called com- 
posite fermiontl. Within mean field theory, where fluctuations of the gauge field 
are ignored, the fermions are assumed to form a conventional Fermi liquid. Many 
authors have considered the stability of the Fermi liquid with respect to fluctua- 
tions of the gauge field. When one calculates the fermionic self-energy E(w) to first 



order in the 
frequencies!*! 



mically screened propagator of the gauge field, one finds for small 



S(w)txw 7 , (1) 



with 7 < 1. This implies a vanishing quasi-particle residue Z — lim w ^o[l — 
9£((jj)/<9aj] -1 , so that the system cannot be a Fermi liquid. Because first order per- 
turbation theory qualitatively changes the analytic structure of the non-interacting 
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Green's function and there is no obvious small parameter in the problem, it is 
necessary to address the effect of higher orders in perturbation theory. 

Clearly, in order to understand the physics of gauge fields in these strongly 
correlated systems, conventional perturbative many-body theory is not sufficient. 
At least one should sum properly chosen infinite sub-classes of Fcynman diagrams. 
Over the past few years a variety of methods have been jropja s e rt to resurn the 
most divergent terms to all orders in perturbation theorju#ElllJOE3llj'oli3llj'0. 
However, so far a general agreement as not been achieved. In particular, in RefO 
it was found by means of a non-perturbative functional integral approacfi00 that 
higher orders completely change the scenario suggested by a one-loop calculation: 
the non-analyticities suggested by lowest order perturbation theory were found to 
be partially removed, so that the spectral function exhibits a well-defined quasi- 
particle peak. This implies that there must exist some higher order terms in the 
perturbative expansion of the self-energy that are more singular than the one-loop 
result. In this work we shall explicitly calculate the two- loop correction to the self- 
energy due to the leading radiative correction to the fermion - gauge field vertex, 
and show that it contains an additional factor proportional to ln(u;o/a>) as com- 
pared with the one- loop result (luq is some finite frequency, see Eq.([72|) below). For 
reasons that will be explained in detail in Sec. 3, this logarithmic correction has been 
missed in a previous calculation by Altshuler, Ioffe and Millisa. The existence of a 
logarithmically divergent correction implies that the true infrared behavior of the 
fcrmionic self-energy cannot be obtained from a one-loop calculation. In fact, the 
prefactor of the logarithm turns out to be negative (relative to the one-loop self- 
energy, see Eq. (f72|) ) , so that it is consistent with the scenario that the summation 
of vertex corrections to infinite orders leads to an_mcrease of the exponent 7 in 
Eq. ([!]), and possibly restores Fermi liquid behaviorEZl (which requires 7 > 1). 

2. The self-energy without vertex corrections 

Let us begin with a careful discussion of the one-loop fcrmionic self-energy correc- 
tion Si due to fluctuations of the gauge field. The relevant Feynman diagram is 
shown in FieJl. Although this diagram has been calculated previously by many 
authors mm let us calculate it once more, using a particular coordinate system 
where all wave- vectors are measured relative to a fixed point on the Fermi surface. 
The insights gained from this calculation will be useful for the more difficult two- 
loop calculation. Using the Matsubara formalism, Fig.[l] represents the following 
expression for the fermionic self-energy, 

= "^E /l f A ^ G o(fc + ?) , (2) 

where V is the volume of the system, (3 is the inverse temperature, and we have 
defined collective labels k = \k,iu> n ], q = [q, iui n ], where Q n = 2n(n + |)//3 are 
fermionic Matsubara frequencies, and to n = 2wn/P are bosonic ones. The super- 
script a indicates that all wave-vectors are measured with respect to a point k a on 
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k+q 



Fig. 1. Leading fluctuation correction to the fermionic self-energy. The solid arrow is the mean- 
field Green's function, and the wavy-line is the RPA screened propagator of the gauge field, see 
Eqs.(S) and (ft. 




Fig. 2. Definition of circular coordinates centered at point k a on the Fermi surface. 



the Fermi surface, as shown in Fig.g. The free Green's function is 

1 



G?(*0 = 7 



(3) 



where £k°+k = v Q ■ k + k 2 /(2m a ) is the mean field energy dispersion (measured 
relative to the chemical potential) for wave- vectors k close to k". Here v° is the 
local Fermi velocity with magnitude Vp , and m a is the local effective mass close to 
k Q . Note that we are not assuming that the Fermi surface is spherically symmetric; 
in particular, for m a — > oo the Fermi surface becomes locally flat. The dynamically 
screened gauge field propagator is within the random phase approximation (RPA) 
and for frequencies in the regime \u> n \ vpq given by 



, RPA, a 



2tt 



[l-(v"-q) 2 ]: 



v F q 



(4) 



where we have used the Coulomb gauge, v° is a unit vector parallel to v° , and the 
energy scale T q is 

T q = v F q {q/qcf = v F q c (q/q c ) 1+n ■ (5) 

The mass m* is some effective mass such that m*vp = kp is a measure for the 
average curvature of the Fermi surface. In contrast, m a vp = k c measures the local 
curvature of the Fermi surface. Throughout this work we shall assume that the 
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momentum scale q c in Eq.(^) is small compared with kp and k c , and that the 
exponent r\ is in the interval 1 < rj < 2. The case rj — 2 is relevant to gauge theories 
of high-temperature superconductors, as well as half- filled quantum Hall systems 
with short-range density-density interactions. 

To perform the integration over the momentum-transfer q in Eq.(|2|), we choose 
the circular coordinates shown in Fig^. For simplicity, we shall restrict ourselves 
to external wave- vectors of the form k = k\\\r a , so that 

k 2 

ek-+k+ q = efc„ + (i + ^ ! -)«F«cosi? + ^, (6) 

where ^ — vpk\\ + k 2 /(2m a ). Keeping in mind that the form ([|) for the gauge 
field propagator is valid for \uj n \ <, v F q, and imposing an ultraviolet cutoff k on the 
g-integration in Eq. (|^) (anticipating that the leading behavior for small frequencies 
is dominated by the infrared singularities of the integrand, we may choose q c <C 
K <C kp) we obtain from Eq.(^) for V — > 00 and /3 — > 00, 

1 f K f VFq v F q 
E"(fc|i,ia)„) = - / dqq / du- 



(2^) 2 m* Jo ™J_ Vrq H+T v 



sin 2 <d 



d$ z =- . ( 7 ) 

■ G^+iu;-(l + ^)v F qcos^-^ 

with Gq 1 = iu n — . Scaling out a factor of Vpq(l + -|^) and defining the complex 
variable 

W M = G ° l+ '"-^ , (8) 
v F q(l + ^) 

Eq.(^) can also be written as 

Ef(*||,iw n ) = j—/ dgg / do; ——— . (9) 

27rm*(i + J) io y-^jpg M + r, 

Here the complex function Z(W) is defined by@ 

= - r * , = w - Vw* - 1 , (10) 

■k J W — cos v 

where the root has to be taken such that \Z\ < 1. It turns out that the leading 
infrared behavior of Eq.(||) is due to the regime where |W(g, u>)\ < 1. To see this, 
note that 

( -isgn{lmW) + W + 0{W 2 ) if|W|<l 
Z(W) = \ . (11) 

[ (2W)- 1 + 0(W- 2 ) if|W|>l 

Because \oj\/(vpq) ;$ 1 and < Vp q m the domain of integration in Eq.(^|), the 
condition \W\ -C 1 is equivalent with 

q > k = max{|fe|||, \u) n \/v F } . (12) 
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To evaluate the integral in Eq.(^), we subdivide the ^-integration into the regimes 
< q < fco and feo < q < k, and use the corresponding asymptotic forms of Z{W) 
given in Eq.(pd|). Using the fact that |i>(W)| < 1, it is not difficult to show that 
the contribution from the regime < q < ko is 

k 2 

s i(fc||.^n)(?<feo) « ^7 ln 0? c / fc o) ■ (13) 

Obviously, for small Q n and fen this correction is negligible. On the other hand, in 
the regime ko < q < k we may approximate (see Eqs.(||) and (pi])) 

Z(W(q, w)) m -isgn(cD„ + u) . (14) 

The corresponding contribution to the self-energy is for |feii | -C k c 

rVFq cL Sm ^ n + Uj ) 

-v F q W+Fq 



E?(fc|| ,«*,)(,>*,) = f—; dqq du . (15) 



Using the fact that the rest of the integrand is an even function of w, we may replace 
under the integral sign 

Sgn(w„ + uj) -> sgn(w„ + u) - sgn(cj) = 2sgn(£„)e(-cj(a>„ + u>)) . (16) 

Then it is easy to show that 

r~ ^ -isgn(d>„) f" ( |u) n |\ 

S 1 (fc||, i w„) (g>feo) = — ^ dqq\n[l + — j . (17) 

Introducing the new integration variable y = T q /\ui n \ and dimensionless frequencies 
[u n — u) n / {vpq c ) and wave-vectors ku = kn/q c , we have 

q = q c \u n y\ 1+ " , = g c w n !+" — — dy , (18) 

1 + 7? 

so that 

Sf(fc||,ia; n )( 9>i . ) . 2 - ._ 

ii - '- = -iff sgn(w„) w„ i+i i?i (fcii, zw n ) , (19) 

where g = q c /kp <C 1 is a dimensionless coupling constant. The dimensionless 
function F\(k\\,iu) n ) is defined by 

Fi (fe|| , iu)„) = 1 / ' dyy ^ ln(l + 1/y) , (20) 
tt(1 +r?) 7 ao 

with the lower limit 

2/o = ^= (^] t^M , (21) 

| fe|| | ^IcDn | - 1 if |o) n |<|fe||| 



u n \ \q c / |W, 
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and the upper limit 



2/i 



1-H7 



(22) 



Note that for large y the integrand in Eq. (p0[) vanishes as y 1 ^ - , so that for 77 > 1 
the integral is ultraviolet convergent. To obtain the limiting behavior of the self- 
energy for \u) n \ — > 0, we may let yi — > 00. Furthermore, for \u> n \ > we may set 
yo = in the lower limit of Eq.(|2C|), so that 



Fi(k\\,iuj n 



dyy ln(l + 1/y) 



Hence, to leading order 

E?(k\\,iu> n )( q>ko ) 



VFq c 



-ic v g sgn(u) n )|u; n | ^ , |w„| > 



(23) 



(24) 



For completeness, let us also discuss the regime |w„| <C where y ^> 1. Then 
F 1 (k h iQ n ) - l^Vo 1 + 0(% 2 ) . (25) 

so that 



S?(fc||,i<D„) (9>feo) 



.9 



tt(77-1) |fc||| 1+ '> 



. |w n |<|*;||| 



(26) 



The spectral function can now be obtained from the retarded self-energy T,f(kn,u>- 
«)+), 



1. 



, w) = Im 



1 



w-&„ -E?(fc||,w + i0+) 



To perform the analytic continuation (i.e. i£D n — > w + i0 + , sgn(tD n ) 
|tD n | — > — zw) we write 



isgn(w„)|<2>„| 1+1 = iw n [-isgn(u)„)iu) n ] J +" 



(27) 
1, and 

(28) 



so that 



isgn(£ n )|w n |i+» 



sgn(u;)|u;| exp 



. .7TJ7-1 
2 7/ + 1 



[A'sgn(w) + zA"] , 



(29) 



with A' = cos(§ 2=1) and A" = sin(f 2=1). Note that for 77 > 1 both A' and A" 
are positive real constants. We conclude that the retarded self-energy is for small 
frequencies in the regime \u>\ > |fcii| given by 



£?(fc||,cj + i0+) 



-c^</|u;| [A'sgn(w) + iA"] 



(30) 
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Note that ImSf(A;|| , u + i0 + ) < 0, so that the retarded Green's function is analytic 
in the entire upper half of the complex frequency plane (as it should), and the 
corresponding spectral function ( p7| ) is positive (as it should). However, the spectral 
function does not exhibit a well-defined quasi-particle peak, because according to 
Eq. (|30|) the imaginary part of the self-energy has the same order of magnitude as 
the real part. Thus, lowest order perturbation theory suggests that fluctuations of 
the gauge field completely destroy the Fermi liquid state predicted by mean field 
theory. 

From our rather detailed derivation it is now easy to identify the regions in 
energy-momentum space that are responsible for this non- Fermi liquid behavior. 
Obviously, in deriving Eq.(|3(]) we may restrict ourselves to the regime where |W(g, w)| 
is small compared with unity, which requires fco q. Moreover, the leading term in 
the expansion of the function Z(W) (see Eq.(|l4[)) for small \ W\ is due to angles •& in 
the integral representation ( |To| ) of this function such that | cos$| \W\ <C 1. This is 
most easily seen by shifting x — — 7r/2in Eq.([To|), and using the fact that for small 
| W\ the leading behavior of the integral is determined by \x — ReW| ^ |ImW| <C 1. 
Hence, 

^^-^I/ x-J-ikr (31) 

Using Im(a; ± ie)^ 1 = =Fi7rsgn(e), we see that Z(W) ~ — isgnlm(VT) for small \W\, 
in agreement with Eq.(nJ). 



The result (30) remains correct if we sum in addition all diagrams where the 
gauge field propagators do not intersect (thus neglecting vertex corrections). This 
is formally achieved by replacing the non-interacting Green's function Gg on the 
right-hand side of Eq.(j|) by the one-loop renormalized Green's function Gf, thus 
turning Eq.(^) into an integral equation. Assuming that the momentum-dependence 
of the so-defined self-consistent self-energy E" is again negligible, the expression of 
E" is still of the form (||), except that W{q,to) should be replaced by 



G n 1 + iu — E?(i<D„ + iu>) — 
W(g, u) = -5 Y \, ■ (32) 



Because sgnlm(PU) = sgnlm(M / ), the resulting self-consistent self-energy is again 
given by Eq.jT^), but with the lower limit ko replaced by fco ~ Qcl&nl 1 ^ > so that 
\W\ £ 1. This leads to the new lower limit yo = \u) n \ in Eq.(p0|). Obviously, for 
\u) n \ — * we may still set yo = 0. Thus, higher order diagrams without crossings of 
the gauge field lines do not modify the result predicted by lowest order perturbation 
theory. 

Let us summarize what we have learned so far. Within a perturbative calcula- 
tion to first order in the RPA screened gauge field propagator, one finds that the 
fluctuations of the transverse gauge field give rise to a non-analytic contribution to 
the fermionic self-energy, implying the non-existence of well defined quasi-particles. 
Denoting by q the momentum transfer and by u> the energy transfer mediated by 
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Fig. 3. Skeleton diagram for the exact self-energy. The solid arrow is the exact Green's function, 
the wavy line is the exact screened gauge field propagator, and the solid triangle represents the 
exact three-legged vertex. 

the gauge field, the non-analyticity in the self-energy arises from energy-momentum 
transfers in the regime 

|cosi?| < M < i ) (33) 
v F q 

where $ is the angle between q and the local Fermi velocity v Q . In other words, the 
non-Fermi liquid behavior is entirely due to momentum transfers that are almost 
parallel to the Fermi surface, corresponding to | cos?9| « |i? ± ?| -C 1- As will be 
discussed in Sec. 3. 2, this greatly reduces the usefulness of the well-known Ward- 
identities for estimating the importance of vertex corrections. 

3. Vertex corrections 
3.1. General remarks 

Given the fact that the leading fluctuation correction qualitatively modifies the 
mean field result, one should worry about higher order corrections. This problem 
has been addressed previously by Altshuler, Ioffe, and MillisEj (AIM), and by Stern 
and HalperinEil. We shall comment on these works below. 

The general structure of the self-energy is conveniently represented in terms of 
the skeleton diagram shown in Fig.[|. The corresponding analytic expression is 

Z a Ck) = -^Y,h?A a Ck;q)G a Ck + q), (34) 

where /i~ is the exact propagator of the gauge field, G a (k) is the exact Green's 
function, and A a (k,q) is the exact three-legged vertex. Obviously, three physically 
different types of corrections can be distinguished. First of all, there are corrections 
to the gauge field propagator beyond the RPA. Because the interaction mediated 
by the gauge field is most si n g ula r for small momentum transfers, the closed loop 
theorem discussed in Refs.E3n3a guarantees that these corrections are small and 
can he safely ignored. This has been explicitly verified at two-loop order by Kim et 
alt3. The second type of corrections consists of diagrams without crossings of gauge 
field lines; these contribute to the exact Green's function G a (k+q) on the right-hand 
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q 




Fig. 4. Leading radiative correction to the fermion - gauge field vertex. The notation is the same 
as in FigM 



side of Eq.(34). At the end of Sec. 2 we have shown that diagrams of this type do not 
modify the leading infrared behavior of the one-loop self-energy. However, it is a 
priori not clear whether this remains truejhr higher order diagrams involving more 
than one loop. In fact, according to Refs.Qia it is essential to take these corrections 
into account by replacing Go — > G\ in internal loops of higher order diagrams. We 
shall come back to this point below. The third type of diagrams are the vertex 
corrections, which are by definition all diagrams contributing to the vertex function 
A a (fc; q). Naively, one could try to expand the vertex function in powers of the RPA 
gauge field propagator, 

oo 

A a (hq) = l + J2 A nCk;q) , (35) 

n=l 

where A" is the sum of all contributions to the three-legged vertex involving n 
powers of /i RPA:Q . However, this expansion might be ill-defined, because successive 
powers might be more and more singular. In this case non-perturbative methods 
are necessary to resum the perturbation series. Here we shall restrict ourselves to 
the more modest task of evaluating the leading vertex correction A™ shown in Fig.||, 
which is explicitly given by 

A?(fc; q) = -j v Y t hf A > a G%Ck + q)G a {k + q + q>) . (36) 

An attempt to evaluate Eq.(|3^) has been made in Refll. However, as will be 
explained in Sec. 3. 3, these authors have missed the dominant contribution. Before 
embarking on a careful evaluation of Eq.([36|), let us explain with the help of Ward- 
identities why vertex corrections can be expected to play an important role in the 
present problem. 



3.2. Ward-identities 

In the limit of vanishing energy- momentum transfer q = [q, iuj] , the three-legged 
vertex A a (fc; q, iuS) can be related to appropriate partial derivatives of the self-energy 
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S Q (fc). In the dynamic limit we first set q = and then let lu — + (seeRefEa). The 
corresponding Ward-identity reads in our notation 

A%Jk) = lim A a (fc; 0, iu) = 1 - ' . (37) 



0)^0 



d{i£j n ) 



On the other hand, in the static limit, we set first ui = and then take the limit 
q — ► 0. This yields 

A« tat (k) = lim iA° (fc; q, 0) = 1 + — -U- . (38) 
q^o o(u_fA:||) 

It has been argued0 that these Ward-identities imply that vertex corrections are 
negligible in the present problem. The argument is based on the observation that 
the gauge field propagator (|J) is singular in the static limit \lu\ -C vpq, so that 
at the first sight it seems that the vertex correction should be estimated from the 
static Ward- identity (^). Because the first order self-energy ( p4| ) is independent 
of fc||, this would imply that the vertex A a (k;q) in Eq.(|3^) can be savely replaced 
by unity. However, this argument is incorrect, because it ignores the fact that the 
non-analyticity of the first order self-energy is entirely due to momentum transfers 
that are essentially parallel to the Fermi surface. To see this more clearly, let us 
recall that for interactions with dominant forward scattering the following more 
general Ward-identity can be derived^, 

[icj-V a ■q}A a (~k;q) = ^ =r- . (39) 

1 G a (k + q/2) G a (k-q/2) 



The dynamic Ward- identity Eq. ( |37| ) can be obtained as a special case of Eq. (|39|) by 
taking the limits u> — + and q — > such that the ratio 

r = v ° -q = ^g cos ^ (40) 



vanishescil Similarly, the static Ward-identity ( |38| ) is obtained by taking these 
limits such that r — > oo. The crucial point is now that according to Eq.(|3^) the 
ratio VFq cosi^/M is small compared with unity for the relevant energy- momentum 
transfers that are responsible for the non-analyticities in the first-order self-energy. 
Hence, the vertex function should be estimated from the dynamic Ward-identity 
([37j) . If we now substitute the perturbative self-energy ( pl| ) into the right-hand side 
of Eq.(|37]), we see that the vertex actually diverges as pnl 1 ^ for u> n — > (recall 
that r\ > 1). It should be kept in mind, however, that the Ward-identity ( |39| ) has 
been derived by linearizing the energy dispersion in the vicinity of the Fermi surface, 
and therefore does contain information about curvature effects. As already pointed 
out in Refs.LrB and discussed in detail below, the curvature of the Fermi surface 



plays a crucial role in the present problem, so that the Ward- identities (|37 39 ) do 



not have much predictive power. In particular, for a spherical Fermi surface the 
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k k+q k+q'+q k+q 1 k 

Fig. 5. Leading vertex correction to the self-energy. 

above Ward-identities cannot be used to obtain qualitative estimates for the order 
of magnitude of vertex corrections. We shall verify the correctness of this statement 
in the next section by explicitly evaluating the leading vertex correction. 



3.3. The leading vertex correction to the self-energy 

The vertex correction A" in Eq. (|36|) gives rise to the following two-loop correction 
to the self-energy, 

1 



= -pvL. h i A "^ ^°o( k + 9) ■ (41) 

The corresponding Feynman diagram is shown in Fig.^. Note that in two dimensions 
the vertex Af(k; q) is a rather complicated function of six variables. From the eval- 
uation of the one-loop self-energy we expect that the regime where the momentum 
transfers q and q' are almost parallel to the Fermi surface will play an important 
role. Below we shall first reproduce the result of AIM§, and then show that these 
authors have missed the dominant contribution. The reason is rather subtle: it 
turns out that the dominant contribution to the self-energy £2 is determined by a 
sub-dominant contribution to the vertex A". 

3.3.1. Exact manipulations 

We begin with the calculation of the function Af(k; q). Setting for simplicity k = 
/c||V Q and introducing the circular coordinates shown in FigJ^, we obtain 

A" (k\\,iu) n ;q, •&,{<*}) = 

* fc+ , Pig' r' d^ W ' U > A \ (42) 

27n; F m*(l + ^)(l+ fc|l+ ^ os V o i-W T i' + M 

with the dimensionless integral I(W, U, A) given by 

1 r* sin 2 #' 

U > A ) = L ^ ^-cos^-coB^-ArintfO ' (43) 

and 
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u 



Go 1 + %J + iuj-{l + ^)v F q cos $ - ^f- 



A = 



in i feii+gcostJN 

v rf( 1 + JL -ir e — ) 

sin$ 



k, 



(45) 
(46) 



The angular integration in Eq.(|4^) can be performed exactly. The result is 

4 




(1 + W')(l + U)[(B - A) 2 + AC 2 ] 

~ 1-VA _ (1 - B)(l - y/B - C 2 /4) 
I- A 



+ 



C 2 



(1-B) 2 + C 2 



B 



(l-B) 2 + C 2 

{B + A){l + B) 
Ay/B - C 2 /A 



where we have defined 
A = 



W - 1 
W' + l 



B 



U-l 
U + l 



c 



2A 

U+l 



(47) 



(48) 



and the roots have to be taken such that Rev^4 > and Re-^/ B — C 2 /4 > 0. To 
make further progress without resorting to numerical methods, the obvious strategy 
is to identify the regimes in the q' — cj'-plane which dominate the integral and then 
find some simplification of the integrand such that the integration can be carried 
out analytically. This is actually not so easy, because we are eventually interested 
in the self-energy (|4~l|), and a priori we cannot exclude the possibility that, due to 
certain symmetry-related cancellations, sub-dominant contributions to the vertex 
A™ are responsible for the dominant infrared behavior of the self-energy . In 
fact, we shall show shortly that this is precisely what happens. 

3.3.2. The contribution from the regime q/q' = 0(1) 

Let us for the moment assume that the leading infrared behavior of the self-energy 
T,2 in Eq.(^l[) is determined by the regime where both parameters \W'\ and \U\ are 
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small compared with unity. At the first sight this assumption seems reasonable, 
because in this case the integrand in Eq.(^) has two poles in the vicinity of the 
real axis, so that the value of the integral is large compared with unity. We now 
show that with this assumption we can reproduce the result of AIMI3. Later we 
shall demonstrate that this assumption is not correct, the infrared behavior of the 
self-energy in Eq.(^Tj) is in fact determined by the regime where \U\ ^> l! 
Obviously, the condition \W'\ 1 is equivalent with (see also Eq.(|l2|)) 



q' £ k = max{|fe|||, \u)„\/v F } , 



q £ qo = max {k , \u\/vf, q 2 /k c } , 



and \U\ £ 1 requires 



where we have used the fact that according to Eq.(|33|) \u>\ ^ vpq\ cosi?|. 
regime it is easy to show from Eq.(^) that to leading order 



(49) 

(50) 
In this 



I(W',U,A) 



[sgn(Im[/) — sgn(lmW') 



1 



1 



U-W 



U-W' + A 



Actually, Eq.(|T|) can be derived in a much simpler way from Eq. 
fact that for small \W\ and \U\ the integral is dominated by r 



(51) 

43) if we use the 
±7r/2. Then we 



may substitute §' — x'±7r/2, replace sin(a;'±7r/2) 
From Eqs.(E3) and (Mq) we find to leading order 



±1, and proceed as in Eq.(tS 



U-W 



iuj — vpq cosd 



Q 

2m a 



VFq 



(52) 



and 



sgn(Im{7) — sgn(ImW // ) = sgn(cj„ + uj + uj r ) — sgn(u) n + u/) . (53) 
Substituting Eqs.(|5l|-|53|) into Eq. (^2|) , we obtain after a straightforward calculation 

Af (fc||,iw„;g,'!?,iw)(g/ > g ) = 



2?rm* 



dq'q' 



sgn(w„) In 11 



sgn(J) n + w) In ( 1 + 
1 



iu — vpq cos •& 



_2L 

2m" 



^Isin^ 



iu — vpq COS $ - 



1- 

2m" 



(54) 



The subscript (q' > qo) indicates that this is the contribution from the regimes 
and Jsol). Note that by definition q > k . 

Let us pause for a moment and compare Eq.(|54|) with the considerations of 
Sec. 3. 2. Because the Ward-identity (|39j) ignores the non-linear terms in the energy 
dispersion, let us take the limit l/m a — > in Eq.(^4|). Comparing the g'-integral 
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with the corresponding integral in the expression for (see Eq.(|l7j)), and taking 
into account that according to Eq.(p3|) |cj| ^, i>f<7| cos and that the integral is 
dominated by frequencies \lu\ ^ |w n |, it is easy so see that for linearized energy 
dispersion A" is indeed proportional to |(D„| 1 +'>. This is in agreement with the 
arguments presented in Sec. 3. 2: in the limit l/m" — > the order of magnitude of 
the vertex correction can be obtained from the dynamic Ward- identity (j37|), and 
not from the static Ward-identity (|38|). 

However, the curvature terms in the denominator of Eq. J54|) cannot be neglected! 
Substituting Eq.(||) into Eq.(|l]), we obtain 



£?(fc||,iw„)(g'><z ) = 



(2ir) 3 (m*) 2 v F J 



dq duj- 



dq'q' 



'III 



sgn(w„) In ( 1 + ) - sgn(w„ + uj) In ( 1 



x 9 

\Cb n 



dQ- 



sin 2 § 



W — cos ■§ 



1 



Wo — cos i? 



+ 



Wo — cos 



f sintf 



(55) 



where is defined in Eq.(||), and Wq — iu;/ {vFq)~q/ (2fc c ). The angular integration 
is the same as in Eq.(fl3"|), and can be done exactly. For \W\ -C 1 (which requires 
q > ko, sec Eq. (|l2|) ) the integral is dominated by i9 ~ ±tt/2, so that we may again 
set sin ps ±1 and extract the leading behavior of the integral as in Eq.([3l|). We 
obtain 



S2(fc||,iai„)(g/ >90 ) 



"^0 



dq'q' In 



'in 



1 - 



27r 2 (m*) 2 

10 



dqq 



fro 



duj- 



r 



g' 



(Go i ) 2 -(^) 2 



(56) 



In the limit l/m a — > the term qq' /m a in the denominator of Eq.(|5"6|) vanishes, so 
that the integral is proportional to Go = [iw n — Cfcii] -1 - Then the dependence on 

2 

uj n can then be scaled out, and we find EfV^i oc \uj n \ 1 +v /(iu) n — This rather 
singular result is drastically modified by the curvature term qq'/m a . This is most 
easily seen by noting that according to Eq.Jlgj) the term qq' jm a scales as |(D rl | T +*' , 
which for r\ > 1 is much larger than Gq 1 . To leading order we obtain from Eq.(|56|) 



££(fe||,io; n )(g'>go) x 



m 
m* 



\lUlr, 



ko 



(57) 



For a circular Fermi surface, where m* = m Q , this agrees with Eq.(8) of AllVli. 
Obviously Eq.(|5^) is negligible compared with the first order self-energy given in 
Eq.©. 

Because according to Eq.([57|) \T^ 2 \ <C AIM argue that the first order self- 
energy correction should actually be included in all internal propagators. In other 
words, the free Green's functions Gq in Eqs.(|36|) and ( [fl| ) should be replaced by 
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the one-loop corrected Green' s function G". If one repeats the above calculation 
with this renormalized propagator, one findsB XS? oc \u> n \ 1+1 . This is much larger 
than Eq.(^), and has the same order of magnitude as X". However, we shall show 
shortly that the dominant infrared behavior of the two-loop diagram shown in Fig.[5] 
is not correctly given by Eq. (^) , but is in fact logarithmically larger than the one- 
loop self-energy X". Consequently, Eq.(|57|) cannot be used to justify the infinite 
resummation of perturbation theory adopted in Ref .□, where all internal propagators 
in higher order diagrams are replaced by one- loop corrected propagators Gf. AIM 
further justify their approach by noting that the reducible diagram corresponding to 
the term X1G0X1 in the first iteration of the Dyson equation (see Fig. 2b of Ren3) 
is larger than the irreducible diagram shown in Fig.|| While such a point of view 
might be valid within the framework of a large-iV expansion, we do not believe 
that for small N the comparison of irreducible and reducible diagrams is physically 
meaningful: because reducible diagrams contain unphysical poles for frequencies 
close to the non-interacting energy dispersion, their contribution to the self-energy 
can always be tuned to be arbitrarily large. 

3.3.3. The contribution from the regime q q' 

We now show that in the limit of vanishing frequency the leading infrared behavior 
of X2 defined in Eq. (]4l]) is not given by Eq. ( |57| ) , but that the two- loop self-energy 
( PH ) in fact exceeds the one-loop self-energy X" by a logarithmically divergent factor, 

XfO||,iij n ) 



-a„hi(l/|w n ) for \uj n \ -> , (58) 



where is a positive numerical constant proportional to m a /m*. The key obser- 
vation is that so far we have tacitly assumed that the leading infrared behavior of 
Eq.(^TJ) is completely determined by the regime where \U\ 1, which means that 
Q 1 ~ <Zo — maxjfco) M/ u f, q 2 /k c }, see Eq. (|50|) . Let us now check the contribution 
from the opposite regime \U\ ^ 1, equivalent with q' £ qo- From Eq.([l7|) (or more 
simply from Eq.([43|)) we see that in this case 

I(W',U,A)^^-L , \U\>1, (59) 

with Z(W) defined in Eq.©. Note that \U\ £ 1 implies \lu\ > v F q', and hence 
g> q' (keeping in mind that v F q >, Because by construction \W'\ <C 1, it is 
clear from Eqs.(^5|) and (|59| ) that in this regime 

W, U, A) « - -^(a> n+ ^)W t q , < qo . (6Q) 

G + vjJ — vf q cos v — 

Substituting this expression into Eq.(^2|) and performing the w'-integration, we see 
that the contribution from the regime q' qo to the vertex A" is 

Kocf, - q ■ \ ~i sgn(d>„) 



-Kin* Gq 1 + iuj - v F qcosd 



2m° 
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x / cWln(l + p) . (01) 



Let us now assume that \tu\ ^ Q(q,'d) 1 where 

2 

^(q,-d) = max{u F gcos#, — -,v F q c \ui n \~} . (62) 

2m a 

The condition \u>\ > vpqd&nl 1 ^ allows us to neglect the term Gq 1 in the de- 
nominator of Eq. (^IJ) and to replace the upper cutoff of the ^'-integral by infinity. 
Introducing again the rescaled frequencies u) n = uj n /(vFq c ) and u> = u>/(vFq c ) (see 
Sec. 2) we obtain in this regime from Eq.(|6l|) 

2 

\ui n I 1+1 

Ai(fe||,i(2) n ;g,'i?,iw)(g/ < g ) ps -ic^g sga(H) n ) — — — , \uj\ £ Cl(q, i?) , (63) 



where the numerical constant c„ is defined in Eq.(|23|). Because by assumption 
g = q c /kp -C 1 and \u>\ ^ |o) n |i+^, the numerical value of the vertex function ( |63| ) 
is certainly small compared with unity. Nevertheless, the contribution from this 
regime dominates the infrared behavior of the self-energy, because Eq. (|63|) is an 
odd function of u, and has therefore a different symmetry than the bare vertex. 
Substituting Eq.@ into Eq.@ and using the fact that for \w\ ^, r2(g,$) we may 
approximate G^{k + q) m (iui)" 1 , we obtain 

Y,%{k\\,iu) n ) {q , <qo) w -T^f— r y dtfJ_M**** 

VFq ^ ^M^kSl A?{k > iQn] q> tf, iuJ ) {q , <qo) . (64) 

-v F q tU(_L q -\- 

Next we note that for sufficiently small q the energy scale T q is negligible compared 
with q 2 /(2m a ), because T q vanishes faster than q 2 . Let us define the wave-vector 
<7i where both energy scales are identical, 

2^ = ^ = VF<lc (|) ^ ' (65) 

i.e. qi = q c {q c /2k c )~^ . It turns out that the leading contribution to Eq. (|64|) is due 
to the regime q < qy. Because by construction |w| > f^g, $) > q 2 /(2m a ), we may 
then neglect the energy T q compared with \u>\ in the integrand. The w-integration 
is now trivial and generates a factor of fl~ 2 (q, Using the expression ( p4] ) for the 
one-loop self-energy Sf , we obtain 

EJ(fc||,i£D n )( g / <go ) = S"(fc||,za) n )i?(j(D„) , (66) 

with 

«c*^) - -TTW-* r ^ r ■ ( 6? ) 

(2tt)^to* 7 ^ 2 (<7,tf) 
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To leading logarithmic order the integration in Eq.(p7|) gives 

4 m a 



where qi is defined by 



R{iu> n ) 



2m a 



v F q c \uj n \ !+" 



Using Eq.(B5h, we obtain for \uj„ \ — > 



(68) 
(69) 

(70) 

a " = ir 2 (l+r))^* ' (71) 

Because of our rather crude method of estimating the integrals, the numerical value 
of a n in Eq.(|7l|) should not be taken too serious, and remains uncertain by a factor 



R{iu n ) ~ — ^ln(l/|w„|) 



where 



of the order of unity. Combining Eqs.(|66|), (|70f) and (|7l|), and taking into account 
that by symmetry the regime q' 3> q gives rise to an equally large contribution to 
the self-energy (this is obvious from the labels in FigJ|) , we finally arrive at Eq. (|5^) . 

Due to the complexity of the integrations, we have not been able to check 
whether the result ( [58| ) is modified if the internal free propagators Gq in Fig.|| 
are replaced by one- loop corrected propagators G". However, as discussed at the 
end of Sec. 3. 3. 2, we do not believe that such an infinite resummation of perturbation 
theory is a sensible procedure (except perhaps within a large- N expansion^). 

4. Conclusions 

In this work we have shown that at the two-loop order the low-frequency behav- 
ior of the self-energy of fermions that are coupled to transverse gauge fields with 
propagator of the form (0) is given by 



S(iw„) = -ic^-%sgn(w n ) 





l+>7 


VFqc 





1 1 ( 

1 - a r . In — — 



(72) 



where c n and a v are positive numerical constants. The logarithmic correction in 
Eq. (^||) is due to the leading radiative correction to the fermion - gauge field vertex, 
and is the main result of this work. It implies that the correct infrared behavior 
of the self-energy of fermions that are coupled to transverse gauge fields cannot be 
obtained from a one-loop calculation. We disagree in this point with the authors of 
Refsfli, who analyzed this problem within a suitably devised 1/iV-expansion. 

Because the leading vertex correction generates an additional logarithm, it seems 
likely that the higher order vertex corrections will give rise to even higher powers 
of logarithms. If we boldly exponentiate the logarithmic correction in Eq.([7^), we 
find 



-ic^-^-sgn(w n ) 
m 



VFq c 



(73) 
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Note that according to Eq. (j7l]) a n is positive and of the order of unity for a spherical 
Fermi surface. It is therefore tempting to speculate that the summation of vertex 
corrections to all orders will indeed push the exponent 7 to a value that is larger than 
predicted by the one-loop result. In fact, the non-perturbative calculation given in 
Rcf suggests that the true exponent 7 is not smaller than unity, implying that 
the infrared fluctuations of the gauge field do not lead to a destruction of the Fermi 
liquid. This would explain why experimentally composite fermions in the half-filled 
Landau level seem to behave as well-defined quasi-particlesEI 
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